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^ ' We re-evaluate the nuclear matrix element for the proton-proton reaction 
^ ■ p + p — s> + + z/e, which is important for stellar-evolution calculations re- 
ly-^ ■ ferring to stars with masses < 1 Mq and for the solar-neutrino problem. We self- 
O , consistently determine the effect of vacuum polarization on the matrix element 
Q>^ '. by first correcting the low-energy scattering data to account for vacuum polar- 
r-| I ization. We then calculate the proton-proton wave function by integrating the 
P-i| Schrodinger equation with vacuum polarization included. We use improved data 
O ' for proton-proton scattering and for the deuteron wave function. We evaluate the 
• uncertainties that are due to experimental error and estimate those that are due to 
ci • theoretical inadequacies. Without vacuum polarization, we find that the square of 
J> . the overlap integral is = 6.96 with an uncertainty of 0.2% due to errors in the 
^ , experimental parameters and an uncertainty of 1% due to lack of knowledge of the 



shape of the nuclear potentials. We estimate the theoretical uncertainty by using 
six different deuteron potentials and five different proton-proton potentials. Vac- 
uum polarization decreases the calculated value by 0.6^q4%. The complete result 
is A^ = 6.92 X (1 ± 0.002^QQQg) where the first uncertainty is due to experimental 
errors and the second uncertainty is due to theoretical uncertainties. Our value 
of A^ is 2% smaller than the value obtained in 1969 by Bahcall and May. The 
improved calculations of the rate of the pp reaction described here increase slightly 
the predicted event rates for the chlorine and the Kamiokande solar-neutrino ex- 
periments. 
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1. INTRODUCTION 



The cross section for the basic reaction that initiates the proton-proton fusion 
chain, p — > + e"*" + fg, is important for calculations of the stellar evolution 
of main-sequence stars of solar mass or less. In particular, the calculated flux of 
the crucial solar neutrinos depends approximately upon the -2.6th power of the 
square of the pp matrix element (Bahcall and Ulrich 1988). The rate for the pp 
reaction was first estimated by Bethe and Critchfield (1938), and then calculated 
with an effective-range approximation by Salpeter (1952). The results of detailed 
numerical calculations were presented by Bahcall and May (1969); various specific 
calculations and corrections have been pubhshed since then (BroUey 1971; Gari and 
Huffman 1972; Gari 1978; Bargholtz 1979; Guessoum 1988; Gould and Guessoum 
1990; and Carlson et al. 1991). 

In this paper we re-calculate the miclear matrix element for the pp reaction us- 
ing improved data and an explicit and consistent treatment of vacuum-polarization 
effects. We also identify and evaluate the various sources of uncertainty. 

The low-energy cross-section factor for the pp reaction can be written (see 
Bahcall et a/., 1982) 

5.(0)^3.80(^) {^)'{^-^)' X 10-MeV-ba.n.. (1) 

The value of Ga/Gv — 1.2573(28) (Hikasa et al. 1992) used here represents a 
weighted average of five precise modern experiments in which the correlation be- 
tween the electron momentum and the neutron spin or the proton recoil spectrum 
are measured. The quantity 5 — O.Ol^Q QQg (Bahcall and Pinsonneault 1992) is the 
fractional correction to the nuclear matrix element due to the exchange of tt and 
p mesons. We do not re-evaluate 5 in this paper (for a discussion of the mesonic 
corrections, see, e.g., Bhn-Stoyle and Papageorgiou 1965, Gari and Huffman 1972, 
Dautry, Rho, and Riska 1976, Bargholtz 1979, and Carlson et al. 1991). We use 
our current result, = 6.92 for the square of the integral of the overlap of the 
proton and deuteron wave functions, which is about 2% smaller than the previous 
best estimate = 7.08. Small effects due to electromagnetic radiative corrections 
are taken into account by using ft values for nuclear beta decays which have not 
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been altered by theoretical corrections for radiative processes (Bahcall and May 
1968). 

We focus on A^, the square of the overlap integral (at zero energy) where 
(Salpeter 1952; Bahcall and May 1969): 

(2 3 \ ^/"^ 
j Ud{r)up^{r)dr. (2) 

The quantity ap is the pp scattering length, and 7 = {2fiEd)^/'^ is the deuteron 
binding wave number (where fi is the proton-neutron reduced mass and is the 
deuteron binding energy). The function Upp{r) is the radial part of the initial pp 
wave function and Udir) is the radial part of the S'-state deuteron wave function. 

In the most systematic previous study, Bahcall and May (1969) carried out 
calculations for a number of different nuclear potentials and wave functions and 
found A2 = 7.08(1 ± 0.025). Since 1969, there have been several published cal- 
culations of (summarized in Table II of Bahcall and Pinsonneault 1992) with 
specific assumptions about the nuclear potentials; the values obtained in all but 
one case ranged from 6.83 to 7.04. Only the calculation by Gould and Giicssoum 
(1990; Guessoum 1988) gives a larger value (A^ = 7.39) than the best estimate of 
Bahcall and May. (Wc discuss problems with the Gould-Guessoum calculation at 
the end of Section 5.) In some of the previous calculations, the effect of vacuum 
polarization was included or partially included, but was not isolated, which makes 
comparison between the different calculations difficult. 

We first calculate A^ — without including vacuum-polarization effects — using 
improved data for pp scattering and for the deuteron wave function. We inves- 
tigate systematically the uncertainties in our result, which are caused in part by 
errors in the experimental measurements of various input parameters and in part 
by the imprecisely known shapes of the pp and np nuclear potentials. The theoret- 
ical uncertainty is smaller than is suggested by a naive comparison of the different 
values for A^ obtained using the nuclear potentials available in the literature. Each 
potential in the literature predicts slightly different values for the measured quan- 
tities. If we constrain each nuclear potential to reproduce the same values for the 
measured parameters, the apparent differences in the calculated values of A^ are 
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reduced. We find A^^ = 6.96. The uncertainty due to experimental errors 

is about 0.2% and that due to theoretical uncertainties is about 1.0%. 

We then include the effect of vacuum polarization in a self-consistent way by 
correcting the nuclear parameters derived from low-energy pp scattering data and 
by including vacuum polarization in the Schrodinger equation used to calculate 
the pp wave function. The first effect raises the value of A^ by 0.41q'^% and the 
second correction decreases A^ by about 1.0%. Both effects must be included in a 
self-consistent treatment. Therefore, vacuum polarization decreases the calculated 
value of A2 by about 0.6i[|i%. The final result is 

A2 = 6.92x(l±0.002lg:0JJ), (3) 

where the first uncertainty is due to experimental errors and the second is due to 
uncertainties in the nuclear potential. 

In Section 2, we discuss the effective-range approximation and use it to estimate 
the uncertainty in A^ due to the errors in the input parameters that are measured. 
In Section 3, we discuss the pp wave functions, and in Section 4, we discuss the 
deuteron wave functions. Our numerical results are presented in Section 5. In 
Section 6, we discuss the effect of including vacuum polarization in our analysis, 
and in Section 7, we summarize our results and discuss briefly implications for 
predicted solar-neutrino rates. 

2. EFFECTIVE-RANGE APPROXIMATION 

Before presenting our numerical results for A^ it is useful to recall the calcu- 
lation of A^ based on the effective-range approximation (Salpeter 1952; Ellis and 
Bahcall 1968; Bahcall and May 1969). 

If Pp and pd are the pp and deuteron effective ranges, Op is the proton scattering 
length, and 7 is the deuteron binding wave number, then 

w2. 



Ae//(0) = A(0)+(^l^)(pd + Pp), (4) 



where 



A(0) = ivO + | 



^i(x) 



X 



(5) 



and N 77^) (1 — 7Pd)] Here, rj^ is the asymptotic ratio of D- to 5'-state 
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deuteron wave functions, and x = (7-R) • 

The low-energy scattering parameters ap and Pp that are used in Eqs. (4) and 
(5) are those determined by fitting low-energy pp phase shifts measured relative 
to a pure Coulomb potential [see Eq. (11) below]. Published values of ap and 
Pp that appear in the literature are often corrected for various electromagnetic 
and/or strong-interaction effects. These corrected values of ap and pp cannot be 
used here. Using ap = —7.8196(26) fm, pp = 2.790(14) fm (Bergervoet et aL, 
1988), Ed = 2.224575(9) MeV, % = 0.0256(4), and pd = 1.759(5) (Brandenburg et 
aL, 1988), we find 

Alff{0) = 6.975 X [1 - 0.29 {pd - 1.759) - 0.10 {pp - 2.790)] . (6) 

We have exhibited in Eq. (6) the dependence of the result on the effective ranges, 
which are the experimental quantities with the largest uncertainties. The uncer- 
tainty in due to the errors in the experimentally determined quantities is 0.2% 
and comes predominantly — but not exclusively — from the imperfect experimental 
knowledge of pd and pp. 

The effective-range approximation yields a value for consistent with what is 
obtained numerically by integrating the Schrodinger equation (~ 6.96; see below 
and Bahcall and May 1969). The analytic result given in Eq. (6) does not assume 
any specific knowledge of the shape of the nuclear potential, reflecting the fact that 
most of the overlap between the proton and the deuteron wave functions occurs at 
radii large compared with the range of the nuclear forces. Since the wave functions 
at asymptotically large radii are accurately determined by the experimentally de- 
termined quantities 7, r]d, Pd, Pp, and ap (Bahcall and May, 1969), A^ is to a large 
extent determined by the experimental quantities and is therefore insensitive to 
the details of the nuclear interaction. This robustness will be quantified further in 
the following sections (see especially the last columns of Table 1 and Table 2). 

Since the effective-range approximation gives an answer in agreement with the 
numerical integrations, we assume that the effect of small changes in the experimen- 
tal parameters on A^ may be estimated well by the effective-range approximation. 
This assumption has been verified quantitatively by comparing changes in the cal- 
culated A^ produced by varying input parameters in simple wave functions. We 
conclude that the experimental uncertainty in A^ is indeed about 0.2%. 
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3. PROTON-PROTON WAVE FUNCTION 



We begin by summarizing the theory of low-energy pp scattering without in- 
cluding vacuum polarization. The s-wave pp radial wave function Upp{r) satisfies 
the radial Schrodinger equation, 



rf2 



U 



u 



-k'^u, 



(7) 



where R = /Me^ = 28.8198 fm, M is the proton mass, V{r) is the nuclear 
potential multiplied by M/Ti?, and k = Mv/2% is the center-of-mass momentum 
where v is the relative velocity. The first boundary condition is that Mpp(O) = 
and the second condition is that Upp{r) approaches a properly normalized distorted 
Coulomb wave at large r. 

Since V{r) has finite range, Upp{r) approaches an asymptotic limit, 0(r), for 
large radii. In this domain, 0(r) satisfies 



dr^ 



1 



-k' 



This equation is solved by 



cl){r)^Co[Go{kr) + cotSFo{kr)], 



(9) 



where Gq and Fq are the irregular and regular Coulomb wave functions, and 5 is a 
phase shift. The Gamow penetration factor is 



^2 

^0 



2777] 



exp(27r77) — 1 ' 



(10) 



where rj — e^/kv, and v is the relative velocity of the two protons. Note that 
the normalization of is such that 0(0) = 1. For a given potential, V{r), the 
phase shifts for each energy, 5{k), are calculated by setting Upp{0) — 0, numerically 
integrating Eq. (7) well beyond the range of V{r), and then matching the solution 
onto Eq. (9). 
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The scattering length and effective range are determined for the measured 
phase shifts in low-energy pp scattering experiments. The scattering length, ap, 
and effective range, pp, are related to the phase shifts by 

C^k cot 6 + ^h{kR) = -— + Ippk'^ + ■■■, (11) 
K dp 1 

where h{x) = (a;^/3)[l + (x^/lO) + ...] for small x. Since the scattering length is 
given by the limit 

^ limC^A;cot5, (12) 



ap fe— »-0 

only the A: ^ solutions to Eqs. (7) and (8) are needed to determine ap. For A; = 0, 
Eq. (8) is solved by modified Bessel functions, and one finds (using the definition 
of Op) that the A; = limit of Eq. (9) is 

cj>{r)=yiir/R)--y2{r/R), (13) 
ap 

where 

yi{x) = 2^/xKl{2^/x), and y2{x) = ^/{x)Il{2^/x). 
Given V{r), the scattering length is then given by evaluating 

1 _ 2 rKo{2^/F/R) + aVrRKi{2^/^) 



ap R aV7^h{2y/rjR) - rh{2-./7jR) ' 



(14) 



at some r large compared with the range of nuclear forces. Eq. (14) converges to ap 
as r is increased. In practice, a radius r > 20 fm yields a precision better than that 
which is determined experimentally. The logarithmic derivative a — ru'pp{r) / Upp{r) 
is determined by solving Eq. (7) numerically with Upp{0) — 0. 

The effective range pp is defined by the relation 

00 

Pp = 2j [(l)\r)-u\r)]dr, (15) 



and can be determined numerically for a given potential V{r). 
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To obtain the pp wave function, we need to assume a form for the nuclear 
potential that yields the measured values of ap and pp, and then integrate the 
Schrodinger equation. We want our estimate of the uncertainties to be conservative; 
therefore, we make minimal assumptions about the nuclear potential. Following 
the example of Bahcall and May (1969), we consider several plausible functional 
forms for V{r). To fit the two experimental quantities, ap and pp, each potential 
is described by two parameters. 

The nuclear interaction is heuristically represented by an attractive potential 
of depth Vq and range b. We try the following forms for the nuclear potential: (i) 
a square well potential, 

( Vq for r < b, 

nr) = \^ . ; (16) 



for r > b; 

(ii) a Gaussian potential, 

V{r) = VoeM-{r/bf]; (17) 

(iii) an exponential potential, 

V{r) ^Voexp{-r/b); (18) 

(iv) a Yukawa potential, 

V{r)^Vo{b/r)eM-b/ry, (19) 

and (v) a repulsive-core (RC) potential, 

V{r) = Vbexp(-r/6) + 300 MeV 6(0.4 fm - r), (20) 

where is a step function. Three of these shapes (square well, Yukawa, and 
exponential) were used by Bahcall and May (1969). The values of Vq and b that 
give ap = —7.8196 fm and pp = 2.790 (see Section II) fm are listed in Table I for 
the five models, and the potentials are plotted in Fig. 1. 
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4. THE DEUTERON WAVE FUNCTION 



The dcuteron wave function is determined by a procedure similar to that used 
to calculate the pp wave functions. To estimate the sensitivity of to details of the 
neutron-proton interaction, we use six wave functions which fit the static deuteron 
parameters {e.g.^ magnetic and electric quadrupole moments), as well as those 
which are more crucial to the calculation of A^: Efi^ and p^^, [see Eqs. (4) and 
(5), and the following discussion]. We are grateful to R. Wiringa for supplying a 
code that calculates the deuteron wave functions. The wave functions are obtained 
from the Argonne vl4 potential (Wiringa et ai, 1984), the Reid soft-core (RSC) 
potential (Reid, 1968), the Urbana vl4 potential (Lagaris and Pandharipande, 
1981), the super-soft-core (SSC) potential (de Tourreil and Sprung, 1973), and the 
Argonne vl8 potential (Wiringa, 1993); in addition, we also calculate A^ with the 
McGee wave function (McGee, 1966), which was used by Bahcall and May. 

The six wave functions are plotted in Fig. 2, and the values of E^, rjd and 
Pd are given in Table 2. Since the potentials used to calculate the wave func- 
tions were constructed by assuming different values for the deuteron properties, 
none of the potentials precisely reproduces the current values of all of the ex- 
perimental quantities of relevance here. Therefore, the uncertainties in the wave 
functions due to uncertainties in the details of the nuclear interaction are not as 
large as suggested by the (already small) differences in the wave functions plotted 
in Fig. 2. For example, the asymptotic behavior of the wave function is given by 
'^d{f) ~^ A'"(27)-^/^ exp(— 7r), where N and 7 are fixed by E^, and rj^- Since 
these primary quantities are determined accurately by the existing experimental 
data [see discussion following Eqs. (4) and (5)], the behavior of the wave func- 
tion at large radii is not as uncertain as suggested by Fig. 2. Most of the overlap 
integral occurs at large radii. The value of A^ that is obtained by numerically 
integrating the overlap of and Upp, when has an incorrect asymptotic behav- 
ior, is inconsistent with current knowledge. In order to obtain an accurate value 
of A^ that uses the current experimental data, we make small corrections with 
the aid of the effective-range approximation. We noted earfier (end of Section II) 
that the effective-range approximation reliably refiects small changes in A^ due to 
small changes in the parameters. Specifically, we evaluate A^ for a given deuteron 
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potential from the relation, 



A2 



Agjj (experimental parameters) 
Ag^^ (model parameters) 



where A^„^ is the value obtained numerically by calculating the overlap integrals 
with the six deuteron wave functions, and the values of Ag// are obtained from 
the effective-range equations, Eqs. (4) and (5). For the potentials listed in Table 
2, the value of C ranges from 0.979 to 0.996. By using several different potentials, 
we probe the sensitivity of A^ to details of the nuclear interaction. 

5. NUMERICAL RESULTS FOR A^ 

In Fig. 3, we plot the integrand Up.p{r)uci{r) as a function of radius for the five 
assumed pp potentials (using the Argonne vl4 deuteron potential). This figure 
shows that drastic changes in the shape of the pp interaction [see Fig 1] result 
in relatively small changes in the value of the integrand, and the difference is 
significant only in the region r < 5 fm. The integrand is insensitive to the shape of 
the pp potential for radii r > 5 fm. In the region where there is a visible difference 
(r < 5 fm), the wave functions are constrained by the measured effective range 
[see Eq. (15)], so those wave functions which are smaller in the region r < 2 fm are 
larger in the region r > 2 fm (and vice versa). About 40% of the integrand comes 
from the region r < 5 fm and only about 2% comes from the region r < 1 fm. 

The values of A^ obtained using the five different pp interactions are listed in 
Table I. (Here we used the Argonne vl4 deuteron potential.) They range from 
A2 = 6.916 (using the square- well potential) to A^ = 6.979 (using; the Yukawa 
potential) to A^ = 6.939 ( using the strong repulsive core). The difference in A^ 
between the pure exponential potential and the potential that includes a strong 
repulsive core is only 0.3%. Thus we assign a total uncertainty of ±0.5% to the 
value of A^ due to uncertainty in the details of the pp interaction. The exponential 
potential yields a central value of A^ = 6.960. 

Next, we determine the spread in the values of A^ caused by different deuteron 
potentials (using the exponential pp potential) . The results are listed in Table II 
and range from A^ = 6.917 (for the RSC potential) to A^ = 6.988 (for the Argonne 
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vl8 potential). We also assign an uncertainty of ±0.5% to the value of due to 
uncertainty in the deuteron wave function. The Argonne vl4 potential yields the 
central value of A^ = 6.960. 

In the column labeled A'^^^ in Table 2, we list the values of A^ obtained by 
naively inserting the various tabulated deuteron wave functions into the overlap 
integral. These values vary from 6.915 to 7.129, which would naively imply an 
uncertainty of about ±1.5%. However, when the values of A^ are corrected for the 
differences in the assumed values of E^, p^, and rjd, the resulting total uncertainty 
is only ±0.5%. The magnitude of the uncertainty in A^ due to uncertainty in the 
shape of the deuteron potential is similar to that due to uncertainty in the shape 
of the pp potential. 

The differences in the various published values for A^ (see Table II in Bahcall 
and Pinsonneault 1992) are primarily the result of differences in the input param- 
eters. The theoretical uncertainty is smaller than one would estimate by simply 
comparing the published results. For example, the result given here is almost 2% 
smaller than that obtained by Bahcall and May (1969). This difference can be 
traced primarily to the high value of embodied in the deuteron wave function 
used therein (from McGee 1966, who adopted E^ = 2.267 MeV), and secondarily 
to smaller differences in the input 7]^, p^, pp, and Qp. If we correct for these dis- 
crepancies in input data using the effective-range approximation, the Bahcall-May 
result becomes 7.02, within 1% of the current result. 

The anomalously high value of A^ = 7.39 obtained by Gould and Guessoum 
(1990; Guessoum 1988) deserves special comment. The authors state that this 
result was obtained using the overlap of the Paris deuteron wave function (Lacombe 
et al. 1980) and a pp wave function employing a Bargmann nuclear potential (see 
Noyes 1967) modified to include a soft repulsive core. Prom the discussion in the 
papers by Gould and Guessoum, it not clear what form they assumed for the 
pp potential, nor is it clear what values of Op and pp they adopted. Therefore, 
it is difficult to make an unambiguous comparison with their result or even to 
understand precisely the origin of the difference between the Gould-Guessoum 
value and all the other values given in Table 1 and Table 2 (as well as the other 
values obtained by different authors; see Table II of Bahcall and Pinsonneault 
1992). We have, however, carried out an illustrative calculation using a potential 
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that seems to be suggested by the description in Gould and Guessoum. We obtain 
a value of 7.42 for A^, consistent with the Gould- Guessoum value, by adding to 
our exponential potential an attractive square-well core of depth 30 MeV and 
range 0.4 fm. However, for this potential, we find values of ap = —8.4826 fm and 
Pp = 2.727 fm, which are inconsistent with the experimental data (discrepancies 
> 200(7 and 4a for ap and pp, respectively). We do not claim that this procedure 
was used by Gould and Guessoum, but it is the only procedure we could think 
of that reproduces their result. In all of the numerical experiments that we have 
performed using wave functions and potentials that are consistent with the current 
experimental data, we have never obtained a value for greater than 7.00. We 
conclude that the Gould- Guessoum value can only be obtained by using some input 
data or some assumption that contradicts the existing experimental information 
on the pp system. 



We next consider, following Bohannon and Heller (1977) and Gould (1991), 
the effect of vacuum polarization (VP) on the pp wave function. In quantum 
electrodynamics, the Coulomb potential is obtained from the Fourier transform 
of the matrix element for scattering via exchange of one virtual photon from an 
electrostatic source. To next order in a, the photon propagator is augmented by 
an electron-positron loop. This augmentation introduces an 0{a) correction to 
the matrix element which, when Fourier transformed, leads to a small correction, 
the Uehling potential (Uehling, 1935), to the electrostatic potential. The complete 
electrostatic potential becomes, in this approximation. 



6. VACUUM POLARIZATION 




(22) 



where 



oo 



1)1/2 



I{r) = J e-2~ (l + ^) 




dx. 



(23) 



1 



The function /(r) has the limiting forms 



7(r) 



7 — 5/6 — In(mer), 



for 



ruer < 1, 



(24) 
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and 

The function /(r) has a logarithmic singularity for very small radii, is of order a few 
until a radius of about l/2me = 193.1 fm, and then suffers an exponential falloff 
(arising from the exchange of a virtual electron-positron pair) for r > l/2me- 

A self-consistent determination of the effect of VP on the calculation of must 
take into account the following two effects. First, VP must be incorporated in the 
analysis of low-energy pp scattering data, which alters the inferred parameters in 
the nuclear potential. Second, the pp wave function must be calculated with the 
Uehling potential in the Schrodinger equation. Using a sophisticated effective- 
range formalism for vacuum polarization, Bohannon and Heller (1977) find that 
including VP in the potential in the Schrodinger equation decreases the pp reaction 
rate by between 0.8% and 1.2%, and using the WKB approximation, Gould (1990) 
finds that VP decreases the pp reaction rate by about 1.3%. However, in both of 
these papers, the quoted result includes only the second effect described above, i.e., 
the inclusion of the Uehling potential in the calculation of the pp wave function. 
Neither paper includes a calculation of the effect of VP on the low-energy scattering 
parameters from which the nuclear potential is determined. 

We perform a numerical calculation of the effect of VP on the pp reaction 
rate which allows us to isolate both effects of VP. First, we consider the effect 
of VP on the low-energy pp scattering parameters, ap and pp, that are obtained 
experimentally. These parameters are fit to phase shifts measured at energies of 
roughly 0.3 MeV to 30 MeV. The effective-range approximation converges only for 
inverse wavenumbers small compared with the range of the nuclear potential. 
Since the range of the Uehling potential is ~ 200 fm, we cannot estimate the effect 
of the Uehling potential on ap and pp by using Eqs. (14) and (15), and incorporating 
the Uehling potential together with the nuclear potential. We must calculate the 
effect of VP on the phase shifts for the energies at which the measurements are 
performed. 

For a given nuclear potential, the phase shifts are determined by integrating 
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the Schrodinger equation, Eq. (7). To account for VP, we make the substitution 

V{r) V{r) 

in Eq. (7). We call the phase shifts obtained from this substitution Svp{k), and 
we call the phase shifts obtained from Eq. (7) {without YP) 6{k). In fact, the mea- 
sured phase shifts are Syp, and the low-energy scattering parameters are correctly 
obtained by fitting to the relation 

Cik cot Svpik) + ^h{kR) = -— + (27) 
K dp 2 

instead of to Eq. (11). The pp potentials used in Section III were constructed 
using scattering parameters obtained from Eq. (11), not Eq. (27). Therefore, in 
constructing the pp nuclear potential as in Section III, we make the substitutions 
Up ^ ttp — Sttp and Pp ^ Pp — 5pp in the low-energy parameters, where Sup and 5pp 
are obtained by fitting the measured phase shifts to 

Ci cot 5vP + ^h(kR) - f(k) = + ^{pp- 5pp)k'. (28) 

K dp — OClp Z 

Here, the function 

f{k) = Clk[cot 5vp{k) - cot 5{k)] (29) 

is obtained for each data point by integrating the Schrodinger equation numerically 
with and without the Uehling potential using a nuclear potential that gives the 
correct Op and pp. We used an exponential nuclear potential, but the results are 
insensitive (to an accuracy of about 0.1% in Op and pp) to the choice of potential. 

We first calculate 5ap and 5pp with the aid of a set of measured phase shifts 
given by Jackson and Blatt (1950). Although the modern data set consists of more 
data points, the Jackson-Blatt data are sufficient (see discussion below of results 
using part of the data and using simulated data) to evaluate the small effect of VP 
on the scattering parameters. By first fitting the complete data set to Eq. (27), 
and then to Eq. (29), we find 5ap = 0.0687 and 5pp = -0.030. Using Eq. (4), 
we find that this results in an increase in of 0.7%. The sign of this result is 
expected: If VP (a repulsive potential) is included, the nuclear potential must be 



2al{r) 
37r 



(26) 
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deeper to compensate, so A is increased. To assess the dependence of our results 
on the data set used, we do the same with only the lower-energy data points. 
(The complete data set consists of 24 data points at ccntcr-of-mass energies of 
0.1765-3.53 MeV obtained with Van dc GraafT generators, and five data points 
with energies of 4.2-14.5 MeV obtained with cyclotrons. In their analysis, Jackson 
and Blatt disregarded the data points obtained with cyclotrons since the data were 
less reliable at the time.) Using only the lower-energy data, we find 5ap = 0.0732 
and 5pp = —0.042, which leads to an 0.4% increase in A^. 

We used simulated data to estimate the effect of VP on the modern data set, 
which include measurements up to 30 MeV. We calculated the effect of VP on the 
values of ap and pp obtained from ten simulated data points that exactly reproduce 
the measured Op and pp in Eq. (11) which are uniformly spaced between 3 and 30 
MeV. This calculation gives the correct fractional changes due to VP if the changes 
are small. For this simulated modern data set, dap = 0.0321 and Spp = —0.001, 
which again results in an increase in A^ of 0.4%. Since this data set most-closely 
resembles the modern data set, we choose 0.4% to be the best estimate of the 
correction to A^ due to VP-corrections to the low-energy scattering parameters. 
To be conservative, we estimate that the total theoretical uncertainty includes the 
full range inferred above using the Jackson-Blatt data and also includes a small 
uncertainty due to the choice of the nuclear potential. We conclude that correcting 
the low-energy scattering parameters for VP results in an 0.4^g^% increase in A^. 

Next we evaluate the effect of including VP on the wave function obtained by 
integrating the Schrodinger equation. By numerically integrating Eq. (7) with and 
without the Uehling potential, but using the same nuclear potential, we find that 
the value of A^ is decreased by 0.9% when the Uehfing potential is included. The 
sensitivity of this result to the choice of the shape of the nuclear potential is less 
than 0.1%. Again, the sign of the effect is expected since the Uehling potential 
is repulsive. This is the result for the correction to A^(0), the matrix element 
squared at zero energy. The most probable energy of interaction in the Sun is 
about 6 keV (Bahcall 1989). We use the WKB approximation to evaluate the 
energy dependence of the VP correction to A^ (Gould 1990; Kamionkowski and 
Bahcall 1993). The magnitude of this VP correction is about 10% larger at 6 
keV than at zero-energy, so we conclude that inclusion of the Uehling potential in 
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the calculation of the wave function decreases by 1.0% at solar energies. Our 
result is consistent with that of Bohannon and Heller, but is slightly smaller than 
Gould's result. This small difference is most likely due to Gould's use of the WKB 
approximation (see Kamionkowski and Bahcall 1993). 

Combining our results from including vacuum polarization in analyzing the 
scattering data and in integrating Schrodinger's equation, we find that VP de- 
creases by 0.6^q'4%. The net effect we find for VP on the pp matrix element 
is about half that found by Bohannon and Heller (1977) and Gould (1991). We 
attribute the difference to the fact that we self-consistently include the effect of 
VP on the measured low-energy scattering parameters, whereas this aspect of the 
influence of VP (which partially cancels the effect of VP in the potential used in 
integrating the Schrodinger equation) was not considered by Gould and was not 
isolated by Bohannon and Heller. 

7. SUMMARY AND DISCUSSION 

We have calculated the matrix element for the reaction p + p ^ ^H + e+ + z/g 
with and without the effects of vacuum polarization. Without vacuum polarization, 
we find the square of the overlap integral to be A^ = 6.96 x (1 ± 0.002 ± 0.010%), 
where the first uncertainty results from la errors in the experimental quantities 
and the second refiects imprecise knowledge of the shape of the nuclear potential. 
We include vacuum polarization self-consistently in the low-energy pp scattering 
parameters and in the numerical calculation of the pp wave function. We find 
that vacuum polarization decreases A^ by 0.6]'^q'4%. Our final result is A^ = 
6.92 X (1 ±0.002^Q'QQg). We show elsewhere that vacuum polarization will decrease 
the rates for all the other reactions in the pp chain and the CNO cycle by small 
amounts, less than 5% (Kamionkowski and Bahcall 1993). 

In addition to the 0{a) VP corrections to the pp potential, there are 0{a) 
radiative corrections (involving an extra soft photon in the final state) to the pp 
reaction. One can use the similarity of the radiative corrections to the axial-vector 
part of neutron decay to those for proton decay in the pp reaction to account 
for these radiative corrections to the pp reaction (Bahcall and May 1968). This 
correction has been included in Eq. (1). 
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It is conventional to use the low-energy cross-section factor, 15(0), in stellar- 
evolution calculations. Following the discussion in Bahcall et al. (1982) and Bahcall 
and Pinsonneault (1992), we find [cf. Eq. (1)]: 

5pp(0) = 3.89 X 10-25 (1 ± 0.011) MeV - barns. (30) 

In calculating the error given in Eq. (30), we have used Icr errors for experimentally- 
measured quantities (such as Ga/Gv or a^) and one-third the total range of values 
for theoretically-calculated effects (such as the uncertainty from nuclear potentials 
or from mesonic corrections). Although the theoretical errors cannot be used in a 
precise statistical sense, the uncertainty quoted in Eq. (30) is intended to be used 
as an approximate la error in Monte Carlo studies {cf. Bahcall and Ulrich 1988) 
of the overall uncertainty in predictions of solar-neutrino event rates. Our best 
estimate is about 3% smaller than that quoted by Bahcall and Pinsonneault. This 
is due in part (2%) to our improved value for A^, and in part (1%) to an updated 
value for Ga/Gv (Hikasa et al. 1992). 

Our results imply a 7.5% increase in the predicted event rates for the Kamiokande 
(Hirata et al. 1991) solar-neutrino experiment and a 6% (0.5 SNU) increase for the 
chlorine solar-neutrino experiment relative to the calculations of Bahcall and Pin- 
sonneault (1992), shghtly increasing the discrepancy between standard-model pre- 
dictions and observations. The predicted event rate for the Borexino experiment 
(sensitive to ^Be, Raghavan 1990) is increased by about 3%. 
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FIGURE CAPTIONS 



1) Proton-Proton potentials. The solid curve is the exponential potential, the 
dot-short-dash curve is the Yukawa potential, the short-dash curve is the 

Gaussian potential, the long-dash curve is the square-well potential, and the 
dot-long-dash curve is the repulsivc-corc potential. All five potentials result 
in a scattering length Op = —7.8196 and an effective range pp = 2.790. 

2) Deuteron wave functions. The solid curve is the wave function for the SSC 
potential, the short-dash-long-dash curve is that for the Urbana vl4 poten- 
tial, the short-dash curve is that for the Argonne vl8 potential, the long-dash 
curve is that for the Argonne vl4 potential, the dot-short-dash curve is that 
for the RSC potential, and the dot-long-dash curve is the McGee wave func- 
tion. 

3) Overlap of pp and deuteron wave functions. As in Fig. 1, the solid curve 
comes from using the exponential pp potential, the dot-short-dash curve from 
the Yukawa potential, the short-dash curve from the Gaussian potential, the 
long-dash curve from the square-well potential, and the dot-long-dash curve 
from the repulsive-core potential. In all five cases, the Argonne vl4 deuteron 
wave function is used. In (a) we show the overlap out to a radius of 50 fm, 
while in (b) we magnify the first 5 fm. 
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Proton-Proton Potentials 


Potential 


Vq (MeV) 


h (fm) 




SW 


-11.751 


2.7718 


6.916 


Gaussian 


-27.729 


1.8912 


6.937 


Exponential 


-98.861 


0.7407 


6.960 


Yukawa 


-46.124 


1.1809 


6.979 


RC 


-314.704 


0.5565 


6.939 



Table 1. Proton-Proton potential parameters. The values of were obtained using 
the Argonne vl4 deuteron potential. 



Deuteron Potentials 


Potential 


Ed (MeV) 


Vd 


Pd (fm) 


A^ 

num 


A^ 


A^ 


Experimental 


2.224575 


0.0256 


1.759 




6.975 




ssc 


2.224066 


0.0255 


1.833 


7.129 


7.128 


6.977 


Urbana vl4 


2.224637 


0.0254 


1.816 


7.099 


7.092 


6.981 


Argonne vl8 


2.224575 


0.0253 


1.770 


7.009 


7.000 


6.988 


Argonne vl4 


2.224884 


0.0266 


1.805 


7.055 


7.070 


6.960 


RSC 


2.224688 


0.0262 


1.758 


6.915 


6.973 


6.917 


McGee 


2.2669104 


0.0269 


1.749 


7.026 


7.031 


6.970 



Table 2. Deuteron potentials. The values of A^ were obtained using the exponential 
proton-proton potential. 
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